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Abstract 

We examine e + e~ — > D~D* + and D*~D+ interactions at 4170 MeV using the CLEO-c detector in 
order to measure the decay constant f D +- We use the Df — > ^ + z/ channel, where the £ + designates 
either a or a r + , when the r + — > 7r + I7. Analyzing both modes independently, we determine 
£(D+ -> fi + u) = (0.594±0.066±0.031)%, and B(Df -» r+z>) = (8.0±1.3±0.4)%. We also analyze 
them simultaneously to find an effective value of B eff (D+ -» = (0.638 ± 0.059 ± 0.033)% and 
extract f D + = 274 ± 13 ± 7 MeV. Combining with our previous determination of B(D + — > /i + ^), 
we also find the ratio f D +/fn+ = 1-23 ± 0.11 ± 0.04. We compare to current theoretical estimates. 
Finally, we find B(D+ e + v) < 1.3 x 10~ 4 at 90% confidence level. 

PACS numbers: 13.20.Fc, 13.66.Bc 
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I. INTRODUCTION 



To extract precise information on the size of Cabibbo-Kobayashi-Maskawa matrix ele- 
ments from B — B mixing measurements, the "decay constants" for Bd and B s mesons or 
their ratio, fB d /fB B , must be well known pj. These factors are related to the overlap of the 
heavy and light quark wave-functions at zero spatial separation. Indeed, the recent measure- 
ment of B® — B° s mixing by CDF j2| that can now be compared to the very well measured B° 
mixing has pointed out the urgent need for precise values [H . Decay constants have been 
calculated theoretically The most promising of these calculations are based on lattice-gauge 
theory that includes light quark loops [||, often called "unquenched." In order to ensure that 
these theories can adequately predict /s s //s d it is useful to check the analogous ratio in 
charm decays f D + / f D + . We have previously measured f D + (l, 0| • Here we present the most 
precise measurement to date of f D + and the ratio f D + / fi>+ . 

In the Standard Model (SM), the only way for a D s meson to decay purely leptonically is 
via annihilation through a virtual W + , as depicted in Fig. [TJ The decay rate is given by [8j 



r 2 

_ f2 



mjM 



VTLc 



M D+ 



\V ri 



(1) 



where M D + is the Df mass, me is the mass of the charged final state lepton, Gf is the Fermi 



coupling constant, and |K< 
we take equal to 0.9738 |3j 



is a Cabibbo-Kobayashi-Maskawa matrix element with a value 




FIG. 1: The decay diagram for 



In this paper we analyze both — > fi + u and — > r + u, r + — > 7r + z7. In both D s 
decays the charged lepton must be produced with the wrong helicity because the D s is a 
spin-0 particle, and the final state consists of a naturally left-handed spin-1/2 neutrino and 
a naturally right-handed spin-1/2 anti-lepton. Because the r + has a mass close to that of 
the Df, the helicity suppression is broken with respect to the fi + decay, but there is an 
additional large phase space suppression. 

New physics can affect the expected widths; any undiscovered charged bosons would 
interfere with the SM W + . These effects may be difficult to ascertain, since they would 
simply change the value of the decay constants. The ratio f D + / is much better predicted 
in the SM than the values individually, so deviations from the the SM expectation are more 
easily seen. Any such discrepancies would point to beyond the SM charged bosons. For 
example, Akeroyd predicts that the presence of a charged Higgs boson would suppress this 
ratio significantly [9|. 
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We can also measure the ratio of decay rates to different leptons, and the SM predictions 
then are fixed only by well-known masses. For example, for t + u to [i + i>\ 



R = 



1 - 



M 



T + 
1 



T(Df -> n+v) 



(2) 



m. 



M 2 , 



Using measured masses [3], this expression yields a value of 9.72 with a negligibly small 
error. Any deviation in R from the value predicted by Eq. [2] would be a manifestation of 
physics beyond the SM. This could occur if any other charged intermediate boson existed 
that affected the decay rate differently than mass-squared. Then the couplings would be 
different for muons and r's. This would be a clear violation of lepton universality [101 ]. 

Previous measurements of f D + have been hampered by a lack of statistical precision, and 
relatively large systematic errors [ll, 12, 13, 14, H, One large systematic error source has 
been the lack of knowledge of the absolute branching fraction of the normalization channel, 
usually D+ — > 07r + [l7j]. The results we report here will not have this limitation. 



II. EXPERIMENTAL METHOD 
A. Selection of D« Candidates 



The CLEO-c detector 18[ is equipped to measure the momenta and directions of charged 



particles, identify them using specific ionization (dE/dx) and Cherenkov light (RICH) 
detect photons and determine their directions and energies. 

In this study we use 314 pb _1 of data produced in e + e~ collisions using the Cornell Elec- 
tron Storage Ring (CESR) and recorded near a center-of-mass energy (Ecu) °f 4.170 GeV. At 
this energy the e + e~ annihilation cross-section into D* + D~+DfD*~ is approximately 1 nb, 
while the cross-section for D+D~ is about a factor of 20 smaller. In addition, D mesons 
are produced mostly as D*D*, with a cross-section of ~5 nb, and also in D*D + DD* final 
states with a cross-section of ~2 nb. The DD cross-section is a relatively small ~0.2 nb 
~20| . There also appears to be DD*tt production. The underlying light quark "continuum" 



background is about 12 nb. The relatively large cross-sections, relatively large branching 
fractions and sufficient luminosities allow us to fully reconstruct one D s as a "tag," and ex- 
amine the properties of the other. In this paper we designate the tag as a D~ and examine 
the leptonic decays of the Df, though, in reality, we use both charges for tags and signals. 
Track requirements, particle identification, 7r°, rj, and K s selection criteria are the same as 
those described in Ref. except that we now require a minimum momentum of 700 MeV/c 
for a track to be identified using the RICH. 

We also use several resonances that decay via the strong interaction. Here we select 
intervals in invariant mass within ±10 MeV of the known mass for rf — > 7t + tt^7], ±10 MeV 
for -> K+K~, ±100 MeV for K*° -> K~tt + , and ±150 MeV for p~ -> t^tt . 

We reconstruct tags from either directly produced D s mesons or those that result from 
the decay of a D*. The beam constrained mass, mec, is formed by using the beam energy 
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to construct the D s candidate mass via the formula 



m BC = y^bcam ~ E ^ ' ( 3 ) 

where % runs over all the final state particles. If we ignore the photon from the D* — > 7_D S 
decay, and reconstruct the mec distribution, we obtain the distribution from Monte Carlo 
simulation shown in Fig. [2j The narrow peak occurs when the reconstructed D s does not 
come from the D* s decay, but is directly produced. 
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FIG. 2: The beam constrained mass mec from Monte Carlo simulation of e + e~ — > DfD*~, — > 
cpTT^ 1 at an Eqm of 4170 MeV. The narrow peak is from the Df and the wider one from D*~ — > jD~ . 
(The distributions are not centered at the or D* + masses, because the reconstructed particles 
are assumed to have the energy of the beam.) 

Rather than selecting events based on only mec, we first select an interval that accepts 
most of the events, 2.015 < wbc < 2.067 GeV, and examine the invariant mass. Distributions 
from data for the 8 tag decay modes we use in this analysis are shown in Fig. [3J Note that 
the resolution in invariant mass is excellent, and the backgrounds not abysmally large, at 
least in these modes. To determine the number of D~ events we fit the invariant mass 
distributions to the sum of two Gaussians centered at the D~ mass, a function we refer to 
as "two-Gaussian." The r.m.s. resolution (a) is defined as 

a = /i<7i + (1 - /i)<7 2 , (4) 

where a± and 02 are the individual widths of each of the two Gaussians and f\ is the fractional 
area of the first Gaussian. The number of tags in each mode is listed in Table [H We will 
later use sidebands of the signal peaks shown in Fig. [3] for part of the background estimate. 

To select our sample of tag events, we require the invariant masses, shown in Fig. [31 
to be within ± 2.5a (±2(7 for the rjp~ mode) of the known D~ mass. Then we look for 
an additional photon candidate in the event that satisfies our shower shape requirement. 
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FIG. 3: Invariant mass of D~ candidates in the decay modes (a) K + K~ir~, (b) K$K~ , (c) r/7r~, 
(d) T]'tt~, (e) 4>p~ , (f) 7r + 7r~7r~, (g) K*~K*°, and (h) rjp~, after requiring the total energy of the 
Dj candidate to be consistent with the beam energy. The curves are fits to two-Gaussian signal 
functions plus a polynomial background. 
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TABLE I: Tagging modes and numbers of signal and background events, within ±2. 5a of the 
D~ mass for all modes, except r/p~ (±2<r), determined from two-Gaussian fits to the invariant 
mass plots, and the number of photon tags in each mode, within ±2.5cr of the D s mass-squared 
determined from fits of the MM* 2 distributions (see text) to a signal Crystal Ball function (see 
text) and a 5th order Chebychev background polynomial, and the associated backgrounds. 



Mode 



Invariant Mass 
Signal Background 



MM* 2 
Signal Background 



K+K--K- 








13871±262 


10850 


8053± 211 


13538 


K S K- 








3122±79 


1609 


1933±88 


2224 


7/7T~; T) — 


' 77 






1609 ± 112 


4666 


1024±97 


3967 


r]'ir-; rj' - 








1196±46 


409 


792±69 


1052 


<t>p~i 4> - 




~K~ 


", p~ — > 7T~7r° 


1678±74 


1898 


1050±113 


3991 










3654±199 


25208 


2300±187 


15723 


K*-K*°; 


K*- 




» K%TT~, K*° - 


* K+tt- 2030±98 


4878 


1298±130 


5672 


vp~; v -» 


77, 


P~ 




4142±281 


20784 


2195±225 


17353 


Sum 








31302 ± 472 


70302 


18645±426 


63520 



Regardless of whether or not the photon forms a D* with the tag, for real D*D S events, the 
missing mass squared, MM* 2 , recoiling against the photon and the D~ tag should peak at 
the mass-squared. We calculate 

MM* 2 = (E CM - E Ds - E,f - {-feu - to. (5) 

where E C m ("P^cm) i s the center-of-mass energy (momentum), E Dg (~P£>J is the energy (mo- 
mentum) of the fully reconstructed D~ tag, and E 1 ( |? 7 ) is the energy (momentum) of the 
additional photon. In performing this calculation we use a kinematic fit that constrains the 
decay products of the D~ to the known D s mass and conserves overall momentum and en- 
ergy. All photon candidates in the event are used, except for those that are decay products 
of the D~ tag candidate. 

The MM* 2 distributions from the selected D~ event sample are shown in Fig. HJ We fit 
these distributions to determine the number of tag events. This procedure is enhanced by 
having information on the shape of the signal function. One possibility is to use the Monte 
Carlo simulation for this purpose, but that would introduce a relatively large systematic 
error. Instead, we use our relatively large sample of fully reconstructed D S D* events, where 
we use the same decay modes listed in Table (TJ we find these events and then examine the 
signal shape in data when one D s is ignored. The MM* 2 distribution from this sample is 



shown in Fig. [5j The signal is fit to a Crystal Ball function 2l|, |22fl . The a parameter, that 
represents the width of the distribution, is found to be 0.032±0.002 GeV 2 . We do expect 
this to vary somewhat depending on the final state, but we do not expect the parameters 
that fix the shape of the tail to change, since they depend mostly on beam radiation and the 
properties of photon detection. 

We fit the MM* 2 distributions for each mode using the Crystal Ball function with fixed 
tail parameters, but allowing a to float, and a 5th order Chebyshev polynomial background. 
We find a total of 18645±426 events within a ±2.5cr interval defined by the fit to each 
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FIG. 4: The MM* 2 distribution from events with a photon in addition to the D s tag for the modes: 
(a) K+K-1T-, (b) K° S K-, (c) rjn- , (d) t/tt", (e) <j>p~ , (f) n+TT'Tr-, (g) K*~K*°, and (h) r\p~ . The 
curves are fits to the Crystal Ball function and a 5th order Chebychev background function. 



mode. There is also a small enhancement of 4.8% on our ability to find tags in (or t + u, 
t + — > n + V) events (tag bias) as compared with generic events, determined by Monte Carlo 
simulation, to which we assign a systematic error of 21% giving a correction of (4.8 ± 1.0)%. 
An overall systematic error of 5% on the number of tags is assigned by changing the fitting 
range, using 4th order and 6th order Chebychev background polynomials, and allowing the 
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FIG. 5: The MM* 2 distribution from a sample of fully reconstructed D s D* + events where one D s 
is ignored. The curve is a fit to the Crystal Ball function. 



parameters of the tail of the fitting function to float. 



B. Signal Reconstruction 

We next describe the search for — > fi + u. Candidate events are selected that contain 
only a single extra track with opposite sign of charge to the tag. The track must make an 
angle >35.9° with respect to the beam line, and in addition we require that there not be any 
neutral cluster detected in the calorimeter with energy greater than 300 MeV. These cuts 
are highly effective in reducing backgrounds. The photon energy cut is especially useful to 
reject — > 7r + 7r°, should this mode be significant, and — > f]'K + . 

Since we are searching for events where there is a single missing neutrino, the missing 
mass squared, MM 2 , evaluated by taking into account the observed fi + , D~, and 7 should 
peak at zero; the MM 2 is computed as 

MM 2 = (Ecu ~ E Da - E 1 - E,f - (V CM - f Ds - f y ~ frf , (6) 

where E^ are the energy (momentum) of the candidate muon track and all variables 

are the same as defined in Eq. 

We also make use of a set of kinematical constraints and fit each event to two hypotheses 
one of which is that the D~ tag is the daughter of a D*~ and the other that the D* s + decays 
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into with the Df subsequently decaying into /i + z/. The kinematical constraints, in the 

center-of-mass frame, are 



In addition, we constrain the invariant mass of the D~ tag to the known D s mass. This 
gives us a total of 7 constraints. The missing neutrino four-vector needs to be determined, 
so we are left with a three-constraint fit. We perform a standard iterative fit minimizing x 2 ■ 
As we do not want to be subject to systematic uncertainties that depend on understanding 
the absolute scale of the errors, we do not make a x 2 cu t but simply choose the photon and 
the decay sequence in each event with the minimum x 2 ■ 

In this analysis, we consider three separate cases: (i) the track deposits < 300 MeV in the 
calorimeter, characteristic of a non-interacting pion or a muon; (ii) the track deposits > 300 
MeV in the calorimeter, characteristic of an interacting pion, and is not consistent with 
being an electron; (iii) the track satisfies our electron selection criteria defined below. Then 
we separately study the MM 2 distributions for these three cases. The separation between 
muons and pions is not complete. Case (i) contains 99% of the muons but also 60% of the 
pions, while case (ii) includes 1% of the muons and 40% of the pions Q. Case (iii) does not 
include any signal but is used later for background estimation. For cases (i) and (ii) we insist 
that the track not be identified as a kaon. For electron identification we require a match 
between the momentum measurement in the tracking system and the energy deposited in 
the Csl calorimeter and we also require that dE/dx and RICH information be consistent 
with expectations for an electron. 

C. The Expected MM 2 Spectrum 

For the fi + i> final state the MM 2 distribution can be modeled as the sum of two Gaussians 
centered at zero (see Eq. HJ). A Monte Carlo simulation of the MM 2 for the <J)tt~ subset of 
K + K~n~ tags is shown in Fig. [6]both before and after the fit. The fit changes the resolution 
from cr=0.032 GeV 2 to cr=0.025 GeV 2 , a 22% improvement, without any loss of events. 

We check the resolution using data. The mode — > K K + provides an excellent testing 
ground. 1 We search for events with at least one additional track identified as a kaon using 
the RICH detector, in addition to a D~ tag. The MM 2 distribution is shown in Fig. [7J 
Fitting this distribution to a two-Gaussian shape gives a MM 2 resolution of 0.025 GeV 2 in 
agreement with Monte Carlo simulation. 

For the t + z/, r + — > tt + V final state a Monte Carlo simulation of the MM 2 spectrum is 
shown in Fig. [SJ The extra missing neutrino results in a smeared distribution. 



1 In this paper the notation K K + refers to the sum of K°K + and K a K + final states. 




(7) 



or E Da 




M Di - M Ds 



143.6 MeV. 



11 



1770107-005 

1000 I i i i i | i i i i | i i i i | i i i i I 




-0.20 -0.10 0.00 0.10 0.20 -0.10 0.00 0.10 0.20 

MM 2 (GeV 2 ) 

FIG. 6: The MM 2 resolution from Monte Carlo simulation for Df — > fj, + u utilizing a (frir~ tag and 
a 7 from either D* decay, both before the kinematic fit (a) and after (b). 

D. MM 2 Spectra in Data 

The MM 2 distributions from data are shown in Fig. [9j The overall signal region we 
consider is —0.05 < MM 2 < 0.20 GeV 2 . The higher limit is imposed to exclude background 
from t]tt + and K°tt + final states. There is a clear peak in Fig. [9](i) due to — > /i + z/. 
Furthermore, the region between the /z + z/ peak and 0.20 GeV 2 has events that we will show 
are dominantly due to the — > t + v decay. The events in Fig. [9](ii) below 0.20 GeV 2 are 
also dominantly due to t + u decay. 

The specific signal regions are defined as follows: for [i + v, —0.05 <MM 2 < 0.05 GeV 2 , 
corresponding to ±2cr; for ru, r + — > n + V, in case (i) 0.05 <MM 2 < 0.20 GeV 2 and in case 
(ii) —0.05 <MM 2 < 0.20 GeV 2 . In these regions we find 92, 31, and 25 events, respectively. 

E. Background Evaluations 

We consider the background arising from two sources: one from real decays and the 
other from the background under the single-tag signal peaks. For the latter, we obtain 
the background from data. We define side-bands of the invariant mass signals shown in 
Fig. [2] starting at 4<r on the low and high sides of the invariant mass peaks for all modes. 
The intervals extend away from the signal peaks by approximately the same width used 
in selecting the signal, 5a, so as to ensure that the number of background events in the 
sidebands accurately reflects the numbers under the signal peaks. Thus the amount of data 
corresponds to twice the number of background events under the signal peaks, except for the 
rjir" and rjp^ modes, where the signal widths are so wide that we chose narrower side-bands 
only equaling the data. We analyze these events in exactly the same manner as those in the 
signal peak. 2 



2 The D s mass used in the fit is chosen to be the middle of the relevant sideband interval. 
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FIG. 7: The MM 2 distribution for events with an identified K + track. The kinematic fit has been 
applied. The curve is a fit to the sum of two Gaussians centered at the square of the K° mass and 
a linear background. 



The backgrounds are given here as the sum of two numbers, the first being the number 
from all modes, except 7]tt~ and i]p~ , and the second being the number from these modes. 
For case (i) we find 2.5+1 background in the p + v signal region and 2.5+0 background in the 
t + u region. For case (ii) we find 2+1 events. Our total sideband background summing over 
all of these cases is 9.0+2.3. The numbers of signal and background events due to false D~ 
tags as evaluated from sidebands are given in Table [TT1 



TABLE II: Numbers of events in the signal region, and background events evaluated from sideband 
regions. 



Case 


Region (GeV 2 ) 


Signal Background 


i 


-0.05<MM 2 < 0.05 


92 


3.5+1.4 


i 


0.05<MM 2 < 0.20 


31 


2.5+1.1 


ii 


-0.05<MM 2 < 0.20 


25 


3.0+1.3 


Sum 


-0.05<MM 2 < 0.20 


148 


9.0±2.3 



This entire procedure was checked by performing the same study on a sample of Monte 
Carlo generated at an Ecm of 4170 MeV that includes known charm and continuum produc- 
tion cross-sections. The Monte Carlo sample corresponds to an integrated luminosity that is 
four times larger than the data. We find the number of background events predicted directly 
by examining the decay generator of the simulation is 28 and the sideband method yields 
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FIG. 8: The MM 2 distribution from Monte-Carlo simulation for Df — > t + v, 
of 4170 MeV. 



~k + v at an Ecm 



22. These are slightly smaller than found in the data, but consistent within errors. We note 
that the Monte Carlo is far from perfect as many branching fractions are unknown and so 
are estimated. 

The background from real decays is studied by identifying each possible source mode 
by mode. For the /i + z/ final state, the only possible background within the signal region is 
Df — > 7r + 7T°. This mode has not been studied previously. We show in Fig. [10] the n + 7i° 
invariant mass spectrum from a 195 pb _1 subsample of our data. We do not see a signal 
and set an upper limit < 1.1 x 10~ 3 at 90% confidence level. Recall that any such events 
are also heavily suppressed by the extra photon energy cut of 300 MeV. There are also some 
Df — > t + u, t + — > ty + V events that occur in the signal region. Using the SM expected ratio 
of decay rates from Eq. [2] we calculate a contribution of 7.4 n + Vv events that we will treat 
as part of the signal. 

For the r + z/, r + — > n + V final state the real Dt backgrounds include, in addition to 



the 7r + 7r° background discussed above, semileptonic decays, possible 7r + 7r°7r° decays, and 
other r + decays. Semileptonic decays involving muons are equal to those involving electrons 
shown in Fig.[9](c). Since no electron events appear in the signal region, the background from 
muons is taken to be zero. The 7r + 7r 7r background is estimated by considering the 7r + 7r + 7r~ 
final state whose measured branching fraction is (1.02±0.12)% [lTJ] . This mode has large 
contributions from /o(980)7r + and other tt + tt~ resonant structures at higher mass 23]. The 
7r + 7r°7r° mode will also have these contributions, but the MM 2 opposite to the 7r + will be at 
large mass. The only component that can potentially cause background is the non-resonant 
component measured by FOCUS as (17±4)% [23]. This background has been evaluated by 
Monte Carlo simulation as have backgrounds from other r + decays, and each is listed in 
Table HD 
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FIG. 9: The MM 2 distributions from data for events with a D s reconstructed in a tag mode, an 
additional positively charged track and no neutral energy clusters above 300 MeV. For case (i) 
when the single track deposits < 300 MeV of energy in the calorimeter. The peak near zero is 
from Df — > [i + v events. Case (ii): track deposits > 300 MeV in the crystal calorimeter but is not 
consistent with being an electron. Case (iii): the track is identified as an electron. 



III. LEPTONIC BRANCHING FRACTIONS 



The sum of MM 2 distributions for case (i) and case (ii), corresponding to the sum of 
— > n + v and — > r + u, r + — > n + V candidates, is compared in Fig. [TT] with the expected 
shape, assuming the SM value of R as given in Eq. [2] for the ratio of t + u to rates. 
The curve is normalized to the total number of events below MM 2 <0.2 GeV 2 . Besides 
the prominent /i + z/ peak and r + u; t + — > tt + V shoulder, there is an enhancement between 
0.25-0.35 GeV 2 , due to K°tt + and r]n + final states, where the decay products other than 
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1770107-007 




FIG. 10: The invariant 7r + 7r° mass. The upper curve shows a fit using a background polynomial 
plus Gaussian signal functions, where the width of the Gaussian is fixed to a value determined by 
Monte Carlo simulation. The lower curve shows just the background polynomial. 



TABLE III: Backgrounds in the Df —>■ t + v, t + — > -k + v sample for correctly reconstructed tags, 
case (i) for 0.05<MM 2 < 0.20 GeV 2 and case (ii) for -0.05<MM 2 < 0.20 GeV 2 . 



Source 




£>(%) # of events case (i) 


# of events case(ii) 


Sum 


Df^ 




8.2 


n+1.8 
u -0 





0+1.8 


Df^ 


vr+7r 7r° 


1.0 


0.03±0.04 


0.08±0.03 


0.11±0.04 


Df^ 




6.4 








T+ 


— > 7r + 7T°I7 


1.5 


0.55±0.22 


0.64+0.24 


1.20±0.33 


T+ 


— > ii^Vv 


1.0 


0.37±0.15 





0.37±0.15 


Sum 






1.0^ 


0.7±0.2 


L - ' -0.4 



the 7r + escape detection. The data are consistent with our expectation that the region 
—0.05 <MM 2 <0.2 GeV 2 contains mostly signal. Recall there are 148 total events only 10.7 
of which we estimate are background, 9.0 from fake D~ tags and 1.7 from real tags and 
decays. Above 0.2 GeV 2 other, larger backgrounds enter. 

The number of real [i + v events N^ u is related to the number of events detected in the 
signal region N dct (92), the estimated background N hkgrd (3.5), the number of tags, N tag , and 
the branching fractions as 



Nuu = Ndct — Vbkgrd — Nt aR • £ 



e'B(D+ -> fji+u) + e"B(D+ -> r+u; t 



7T V , 



where e (80.1%) includes the efficiency for reconstructing the single charged track including 
final state radiation (77.8%), the (98.3±0.2)% efficiency of not having another unmatched 
cluster in the event with energy greater than 300 MeV, and for the fact that it is easier 
to find tags in events than in generic decays by 4.8%, as determined by Monte Carlo 
simulation. The efficiency labeled e' (91.4%) is the product of the 99.0% muon efficiency 
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FIG. 11: The sum of case (i) and case (ii) MM 2 distributions (histogram) compared to the predicted 
shape (curve) for the sum of — ► /i + z/ and r + z^, r + — » tt + F. The curve is normalized to the total 
number of events below MM 2 <0.2 GeV 2 . 



for depositing less than 300 MeV in the calorimeter and 92.3% acceptance of the MM 2 cut 
of |MM 2 | < 0.05 GeV 2 . The quantity e" (7.9%) is the fraction of r+z/; r+ -> vr+17 events 
contained in the fi + u signal window (13.2%) times the 60% acceptance for a pion to deposit 
less than 300 MeV in the calorimeter. 

The two branching fractions in Eq. [8] are related as 

B{D+ -> t + u; t + -> 7T+17) = R-B(t + -> ti + V)B{D+ -> //+!/) = 1.059 -> ^ + z/) , (9) 

where we take the Standard Model ratio for i? as given in Eq. [2] and *B(t + — > 
7r +I7) = (10.90±0.07)% 0- This allows us to solve Eq. [HJ Since iV det = 92, A^ bkgrd =3.5±1.4, 
and A^ ta g = 18645 ± 426 ± 1081, we find 

B{Df -> n + v) = (0.594 ± 0.066 ± 0.031)%. (10) 

We can also sum the [i + v and r + v contributions, where we restrict ourselves to the MM 2 
region below 0.20 GeV 2 and above -0.05 GeV 2 . Eq. [8] still applies. The number of events 
in the signal region and the number of background events changes to 148 and 10.7jr|; 3 , 
respectively. The efficiency e' becomes 96.2%, and e" increases to 45.2%. Using this method, 
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we find an effective branching fraction of 

B eS {D + -> fj, + u) = (0.638 ± 0.059 ± 0.033)%. (11) 

The systematic errors on these branching fractions are given in Table IIVI The error on 
track finding is determined from a detailed comparison of the simulation with double tag 
events where one track is ignored. "Minimum ionization" indicates the error due to the 
requirement that the charged track deposit no more than 300 MeV in the calorimeter; it 
is determined using two-body D° — ► K~n + decays (see Ref. 0]). The error on the photon 
veto efficiency, due to the 300 MeV/c extra shower energy cut, is determined using Monte 
Carlo simulation. The Monte Carlo was cross-checked using a sample of fully reconstructed 
D + D~ events and comparing the inefficiency due to additional photons with energy above 
300 MeV/c. These events have no real extra photons above 300 MeV/c; those that are 
present are due to interactions of the D+ decay products in the detector material. The error 
on the number of tags of ±5% has been discussed earlier. In addition there is a small error 
of ±0.6% on the t + u branching fraction due to the uncertainty on the r + decay fraction 
to 7i + V. Additional systematic errors arising from the background estimates are negligible. 
Note that the minimum ionization error does not apply to the summed branching fraction 
given in Eq. [TTJ in this case the total systematic error is 5.1%. 



TABLE IV: Systematic errors on determination of the D + — > [i+v branching fraction. 



Error Source 


Size (%) 


Track rinding 


0.7 


Photon veto 


1 


Minimum ionization* 


1 


Number of tags 


5 


Total 


5.2 



*-Not applicable for summed rate 

We also analyze the r + v final state independently. We use different MM 2 regions for 
cases (i) and (ii) defined above. For case (i) we define the signal region to be the interval 
0.05<MM 2 <0.20 GeV 2 , while for case (ii) we define the signal region to be the interval 
-0.05<MM 2 <0.20 GeV 2 . Case (i) includes the /i + z/ signal, so we must exclude the region 
close to zero MM 2 , while for case (ii) we are specifically selecting pions so the signal region 
can be larger. The upper limit on MM 2 is chosen to avoid background from the tail of the 
K°7r + peak. The fractions of the MM 2 range accepted are 32% and 45% for case (i) and (ii), 
respectively. 

We find 31 signal and 3.51};! background events for case (i) and 25 signal and 5.1±1.6 
background events for case (ii). The branching fraction, averaging the two cases is 

B(D+ -> t + u) = (8.0 ± 1.3 ± 0.4)%. (12) 

Lepton universality in the Standard Model requires that the ratio R from Eq. [2] be equal to 
a value of 9.72. We measure 

fl ' rkw^r 13 - 4±26±0 - 2 - (13) 
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Here the systematic error is dominated by the uncertainty on the minimum ionization cut 
that we use to separate the /z + z/ and t + u regions at 300 MeV. We take this error as 2%, since 
a change here affects both the numerator and denominator. The ratio is consistent with the 
Standard Model prediction. Current results on D + leptonic decays also show no deviations 
26[. The absence of any detected electrons opposite to our tags allows us to set an upper 



limit of 

B(D+ -> e + u) < 1.3 x 10~ 4 (14) 

at 90% confidence level; this is also consistent with Standard Model predictions and lepton 
universality. 

IV. CHECKS OF THE METHOD 

We perform an overall check of our procedures by measuring B(D+ — > ~K°K + ). For this 
measurement we compute the MM 2 (Eq. [H]) using events with an additional charged track 
but here identified as a kaon. These track candidates have momenta of approximately 1 
GeV/c; here our RICH detector has a pion to kaon fake rate of 1.1% with a kaon detection 
efficiency of 88.5% [19]. For this study, we do not veto events with extra charged tracks, or 
neutral energy deposits >300 MeV, because of the presence of the K°. 

The MM 2 distribution is shown in Fig. [7J The peak near 0.25 GeV 2 is due to the decay 
mode of interest. We fit this to a linear background from 0.02-0.50 GeV 2 plus a two-Gaussian 
signal function. The fit yields 375±23±18 events. Events from the r]7T + mode where the 
7r + fakes a K + are very rare and would not peak at the proper MM 2 . Since i]K + could in 
principle contribute a background in this region, we searched for this final state in a 195 pb _1 
subsample of the data. Not finding any signal, we set an upper limit of B(D+ — > r]K + ) < 
2.8 x 10~ 3 at 90% confidence level, approximately a factor of ten below our measurement. 
This final state would peak at a MM 2 of 0.30 GeV 2 and would cause an asymmetric tail on 
the high side of the peak. Since we see no evidence for an asymmetry in the K K + peak 
we ignore the r]K + final state from here on. In order to compute the branching fraction we 
must include the efficiency of detecting the kaon track 76.2%, including radiation [24], the 



particle identification efficiency of 88.5%, and take into account that it is easier to detect 
tags in events containing a K K + decay than in the average D S D* event due to the track 
and photon multiplicities, which gives a 3% correction. 3 These rates are estimated by using 
Monte Carlo simulation. We determine 



B(DJ -> ICK+) = (2.90 ± 0.19 ± 0.18)%, (15) 

where the systematic errors are listed in Table [V] We estimate the error from the signal 
shape by taking the change in the number of events when varying the signal width of the 
two-Gaussian function by ±lcr. The error on the background shapes is given by varying the 
shape of the background fit. The error on the particle identification efficiency is measured 
using two-body D° decays [HI . The other errors are the same as described in Table IIVI 
Again, the largest component of the systematic error arises from the number of tag events 
(5%). In fact, to use this result as a check on our procedures, we need only consider the 



The tag bias is less here than in the n + v case because of the K° decays and interactions in the detector. 
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systematic errors that are different here than in the y^v case. Those are due only to the 
signal and background shapes and the particle identification cut. Those systematic errors 
amount to 3.7% or ±0.11 in the branching fraction. 

To determine absolute branching fractions of charm mesons, CLEO-c uses a method 
where both particles are fully reconstructed (so called "double tags") and the rates are 
normalized using events where only one particle is fully reconstructed. Our preliminary 
result using this method for B(L>t — > i^5i^ + ) = (1.50±0.09 ± 0.05)%, which when doubled 
becomes (3.00±0.19±0.10)% [25J. This is in excellent agreement with the number in Eq. [TBI 
These results are not independent. 



TABLE V: Systematic errors on determination of the Df — > K K + branching fraction. 



Error Source 


Size (%) 


Signal shape 


3 


Background shape 


2 


Track rinding 


0.7 


PID cut 


1.0 


Number of tags 


5 


Total 


6.3 



We also performed the entire analysis on a Monte Carlo sample that corresponds to 
an integrated luminosity four 4 times larger than the data sample. The input branching 
fraction in the Monte Carlo is 0.5% for fi + u and 6.57% for t + u, while our analysis measured 
B(D+ -> fj+v) =(0.514±0.027)% for the case (i) fi+u signal and (0.521±0.024)% for (j+u 
and t + v combined. We also find (6.6±0.6)% for the r + v rate. 

V. THE DECAY CONSTANT 

Using our most precise value for B(Df — > /i + z/) from Eq. [TTJ that is derived using both 
our fi + u and t + v samples, and Eq. [I] with a D s lifetime of (500±7) x 10~ 15 s [3], we extract 

f D + = 274 ± 13 ± 7 MeV. (16) 

We combine with our previous result ^ 

f D + = 222.6 ± 16.7±|;| MeV (17) 

and find a value for 

£ 

= 1.23 ±0.11 ±0.04, (18) 

JD+ 

where only a small part of the systematic error cancels in the ratio of our two measurements. 
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VI. CONCLUSIONS 



Theoretical models that predict f D + and the ratio j^j- are listed in Table |VH Our result 
for fo s is slightly higher than most theoretical expectations. We are consistent with lattice 
gauge theory, and most other models, for the ratio of decay constants. There is no evidence 
at this level of precision for any suppression in the ratio due to the presence of a virtual 
charged Higgs |9J|. 



TABLE VI: Theoretical predictions of f D +, and f D +/fo+- QL indicates quenched lattice 



calculations. 



Model 


f D + (MeV) 


f D+ (MeV) 


/d+//d+ 


Lattice (HPQCD+UKQCD) [27] 


241 ±3 


208 ±4 


1.162 ±0.009 


Lattice (FNAL+MILC+HPQCD) [28] 


249 ± 3 ± 16 


201 ± 3 ± 17 


1.24 ±0.01 ±0.07 


QL (QCDSF) [29] 


220 ± 6 ± 5 ± 11 206 ± 6 ± 3 ± 22 1.07 ± 0.02 ± 0.02 


QL (Taiwan) [30J 


266 ± 10 ± 18 


235 ± 8 ± 14 


1.13 ±0.03 ±0.05 


QL (UKQCD) [31] 


236 ± 8±\l 


210 ± 10^ 


1.13 ±0.021°;°:* 


QL [32] 


231 ± 12lf 


211 ± 14±f 2 


1.10 ±0.02 


QCD Sum Rules [33] 


205 ± 22 


177 ±21 


1.16 ±0.01 ±0.03 


QCD Sum Rules [34] 


235 ± 24 


203 ± 20 


1.15 ±0.04 


Field Correlators [35] 


210 ±10 


260 ± 10 


1.24 ±0.03 


Quark Model [36] 


268 


234 


1.15 


Quark Model [37J 


248±27 


230±25 


1.08±0.01 


LFQM (Linear) [38] 


211 


248 


1.18 


LFQM (HO) [38] 


194 


233 


1.20 


LF-QCD [39] 


253 


241 


1.05 


Potential Model [40J 


241 


238 


1.01 


Isospin Splittings [41] 




262 ± 29 





By using a theoretical prediction for f D + / f D + we can derive a value for the ratio of CKM 
elements |V^/V^|. Taking the value from Ref. [13] of 1.162 ± 0.009, we find 

\V cd /V cs \ = 0.2171 ± 0.021 ± 0.0017 , (19) 

where the first error is due the statistical and systematic errors of the experiment and the 
second is due to the stated error on the theoretical prediction. This value is expected to be 
almost equal to the ratio of the CKM elements IK^I/IV^I 

We now compare with previous measurements. The branching fractions, modes, and 
derived values of f D + are listed in Table IVTTl Our values are shown first. We are generally 
consistent with previous measurements, although ours are more precise. 

Most measurements of Df — > t+v are normalized with respect to B(D+ — > (pn + ) = B^. 
An exception is the OPAL measurement which is normalized to the D s fraction in Z° events 
that is derived from an overall fit to heavy flavor data at LEP [42]. It still, however, relies on 
absolute branching fractions that are hidden by this procedure, and the estimated error on 
the normalization is somewhat smaller than that indicated by the error on available at 
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TABLE VII: Our results for B(Df — ► /i + ^), B(Df — > t + z^), and / D + compared with previous 
measurements. Results have been updated for the new value of the D s lifetime 3. ALEPH 
combines both measurements to derive a value for the decay constant. 



Exp. 


Mode 


B 


Bfr (%) 


/„+ (MeV) 


CLEO-c 




(5.94 ±0.66 ±0.31) • 10~ 3 




264 ± 15 ± 7 


CLEO-c 




(8.0 ± 1.3 ±0.4) • 10~ 2 




310 ± 25 ± 8 


CLEO-c 


combined 






274 ± 13 ± 7 


CLEO [11] 




(6.2 ± 0.8 ± 1.3 ± 1.6) • 10~ 3 


3.6±0.9 


273 ± 19 ± 27 ± 33 


BEATRICE [12] y+v 


(8.3 ±2.3 ±0.6 ±2.1) • 10~ 3 


3.6±0.9 


312 ± 43 ± 12 ±39 


ALEPH [13] 


fjL + V 


(6.8 ± 1.1 ± 1.8) ■ 10~ 3 


3.6±0.9 


282 ± 19 ± 40 


ALEPH [13] 


T + U 


(5.8 ±0.8 ±1.8) • 10- 2 






L3 [14] 


T + U 


(7.4 ± 2.8 ± 1.6 ± 1.8) • 10~ 2 




299 ± 57 ± 32 ± 37 


OPAL [15] 


T + V 


(7.0 ±2.1 ±2.0) • 10~ 2 




283 ± 44 ± 41 


BaBar [16] 


y + v 


(6.74 ± 0.83 ± 0.26 ± 0.66) • 10" 


' 3 4.71±0.46 


283±17±7±14 



the time of their publication. The L3 measurement is normalized taking the fraction of D s 
mesons produced in c quark fragmentation as 0.11±0.02, and the ratio of D*/D s production 
of 0.65±0.10. The ALEPH results use B^ n for their n + u results and a similar procedure as 
OPAL for their t + u results. We note that the recent BaBar result uses a larger B^ than 
the other results. 

The CLEO-c determination of f D + is the most accurate to date. It also does not rely 
on the independent determination of any normalization mode. (We note that a preliminary 
CLEO-c result using — > r + z/, t + — ► e + T7u [43] is consistent with these results.) 
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